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ABSTRACT 

- Modelling method for optimizing faster and more realistically the displacement 
conditions, in a porous medium wettable by a first fluid (water for example), of a 
mixture of fluids including this wetting fluid, another, non-wetting fluid (oil for 
example) and a gas. 

- The method comprises experimental determination of the variation curve of the 
capillary pressure in the pores as a function of the saturation in the hquid phases, 
modelling the pores of the porous medium by means of a distribution of capillaries with 
a fractal distribution by considering, in the case of a three-phase water (wetting fluid)- 
oil-gas mixture for example, a stratification of the constituents in the pores, with the 
water in contact with the walls, the gas in the centre and the oil forming an intercalary 
layer, determination, from this capillary pressure curve, of the fractal dimension values 
corresponding to a series of given values of the saturation in the liquid phase, modelling 
the hysteresis effects that modify the mobile saturations of the fluids effectively 
displaced in the sample, that vary during drainage and imbibition cycles, modelling the 
relative permeabilities directly in the form of analytic expressions depending on the 
various fractal dimension values obtained and in accordance with the stratified 
distribution of the various fluids in the pores, and using a porous medium simulator in 
order to determine, from the relative permeabilities, the optimum displacement 
conditions of the fluids in the porous medium. 

- Application : petroleum production, soil depollution, etc. 
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FIELD OF THE INVENTION 

The present invention relates to a method for modelling two-phase or three-phase 
flows in a porous medium, in drainage and imbibition. It is based on a fractal 
representation of the porous medium and on an original approach for handling 
phenomena linked with hysteresis (change in the direction of variation of the saturations). 

The method according to the invention is applicable in many fields where fluid flows 
in porous media are to be modelled in order to optimize the conditions of the 
displacement thereof in drainage and imbibition. Examples of fields of application are : 

a) development of an oil reservoir and notably enhanced hydrocarbon production by 
injection of fluids, using for example alternate iiyections of liquid and gas slugs (a 
method referred to as WAG). It constitutes an advantageous tool allowing reservoir 
engineers to study also well productivity and injectivity problems ; 

b) soil depollution and notably depollution of industrial sites by injection of 
substances such as surfectants in polluted layers; 

c) cleaning of filter cakes by displacement of the substances retained therein; 

d) wood drying ; 

e) optimization of chemical reactions for example by displacement of reaction 
products in a catalyst mass in order to increase the surfaces of contact, etc. 

The method according to the invention is directly applicable by reservoir engineers in 
order to determine, for example, the most suitable enhanced recovery method to be 
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applied to an underground hydrocarbon reservoir. The method can also serve within the 
scope of industrial site depollution operations for example. 

BACKGROUND OF THE INVENTION 
1) Experimental studies 

5 Experimental determination of the relative permeabilities of a porous medium 

wherein a multiphase fluid flows is not easy. Measuring operations are usually simplified 
by considering that one of the phases is immobile In a state of irreducible saturation. 

The values are for example acquired by means of a well-known e>q)erimental method 
referred to as "steady state" for determining relative permeabilities, which consists in 

10 allowing a three-phase fluid to flow with imposed flow rates between the phases. The 
relative permeabilities expressed as a fimction of the two saturations are calculated by 
applying Darcy's law to each phase. It is not an established fact that the relative 
permeability measurements obt^ed by means of this method are really representative of 
the fluid displacements and, in any case, they take a long time because, at each regime 

15 change, one has to wait for a state of equilibrium. 

Another known method consists in carrying out laboratory tests in order to 
determine measurement tables (as shown in Fig.l) relating the relative permeabilities and 
the saturations for each pair of fluids of the three-phase mixture. By adjusting 
experimental production curves, one tries to progressively adjust the three-phase relative 
20 permeabilities. These data tables are then entered iato an Athos® type numerical 
simulator which computes the fluid productions. This method being based on the prior 
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acquisition of many experimental measurements progressively adjusted by calibration, it 
takes a long time. 

2) Relative permeability models 

The known empirical model referred to as Stone's model allowrs, by en^irical 
correlations, to predict data relative to a three-phase flow from data corresponding to a 
two-phase flow. It is valid only in case of a high water wettability and it is generally 
considered to be poorly predictive. 

There are two known types of physical models for modelling three-phase flows, 
based on capillary pressure curves. The capillary pressure curves are connected with a 
saturation (for exanple that of the mercury injected) and a pore radius from which the 

. . _ 2cr 

mercury stops for a given injection pressure, determined by Laplace's law, Pmj. — . 

A first porous media representation model is described by : 

- Burdine, N.T. : "Relative permeability calculations from pore size distribution 
data". Trans AIME (1953), VoL198, or by 

- Corey, A.T. : "The interrelation between oil and gas relative permeabilities". Prod. 
Monthly (1954), VoL19, 38. 

According to this model, the porous medium is represented by a bundle of cylindrical 
capillaries with a radius distribution given by the capillary pressure curve obtained by 
mercury injection. The permeabilities are obtained by applying Poiseuille's law to the 
flow of fluids in these capillaries. 
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This model is based on the representation of the porous medium as an assembly of 
capillaries with different radii. The relation between the volume and the radius of the 
pores is given by the value of the slope of the pseudo-plateau. The three fluids are 
supposed to share the capillaries between them, the wetting fluid (water) occupjong the 

5 smallest ones, the least wetting fluid (gas) the largest ones and the third fluid (oU) a zone 
with intermediate-size pores. It is not possible to describe the interactions between the 
fluids because, in such a model, they flow through separate channels. Finally, this model 
can be useful only if the pseudo-plateau covers a wide range of saturations. According to 
this model, the three phases of a three-phase flow move in different capillaries and there 

10 is no interaction between them. 

Another known physical porous medium representation model is described 

by: 

- de Gennes, P.O. : "Partial Filling of a Fractal Structure by a Wetting Fluid", 
Physics of Disordered Materials 227-241, New York Plenum Pub. Corp. (1985), taken 

15 up by 

- Lenormand, R. : "Gravity Assisted Inert Gas Injection : Micromodel Experiments 
and Model based on Fractal Roughness", The European Oil and Gas Conference 
Altavilla MiUca, Palermo, Sicily (1990). 

AccordiQg to this model, one considers that the inner surface of the pores is isotropic 
20 and has a fractal character, and it can be modelled as a "bunch" of parallel capillary 
grooves so that the pores exhibit a fractal cross section. The cross section of each pore is 
constructed according to an iterative process (Fig.l). The half-perimeter of a circle of 
radius Ro is divided into t] parts and each of these rj parts is replaced by a semi-circle or 
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groove. At each stage k of the process, new semi-circular grooves of radius and 
of total section Aj^ are created. 

The fractal dimension DL of the cross section at the end of stage k is related to the 
number of objects Nk generated with the given scale Ij^ by the relation : 



The fractal dimension can be deduced from a mercury capillary pressure curve 
according to the following procedure. Mercury is injected into a porous medium with an 
injection pressure that increases in stages. Laplace's law allows to deduce the pore 
volume, knowing the volume of mercury injected for a given capillary pressure and one 
can construct the drainage capillary pressure ciirve relating the injection pressure to the 
amount of mercury injected and the curve relating the proportion of the total volume 
occupied by the pores and the size of the pores. In cases where a wetting liquid is 
drained from the porous medium such as water by gas injection, the correlation between 
the gas-water capillary pressure and the saturation of the wetting phase is given by : 



The experimental results readily show that the values of the gas-water relative 
permeabilities ospressed as a flinction of the three saturations, obteuned from the 
expressions given by the known models and the phase distribution modes in the structure 
of the pores, are far from the measured values and therefore that the models concerned 
prove to be too simplistic to represent the complex interactions that take place between 
the fluid phases. 
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Pc = S, 
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Patent FR-2,772,483 (US-6,02 1,662) describes a modelling method for optimizing 
fester and more realistically the flow conditions, in a porous medium wettable by a &st 
fluid (water for example), of a mixture of fluids including this wetting fluid and at least 
another fluid (oil and possibly gas). This method involves modelling the pores of the 

5 porous mediimi by a distribution of capillaries with a ifractal distribution considering, in 
the case of a three-phase water (wetting fluid)-oil-gas mixture for example, a 
stratification of the constituents in the pores, with the water in contact with the walls, the 
gas in the centre and the oil forming an intercalary layer. It con^rises experimental 
determination of the variation curve of the capillary pressure in the pores as a fimction of 

10 the saturation in the liquid phases, fi"om which the fractal dimension values 
corresponding to a series of given values of the saturation in the liquid phase are 
deduced. It also comprises modelling the relative permeabilities directly in form of 
analytic expressions depending on the various firactal dimension values obtained and in 
accordance with the stratified distribution of the different fluids in the pores. A porous 

15 medium simulator is used fi"om these relative permeabilities to determine the optimum 
conditions of displacement of the fluids in the porous medium. 

The hysteresis phenomenon relates to the variations in the petrophysical properties 
(relative permeabilities, capillary pressure, resistivity index, etc.) observed according to 
whether measurements are performed during drainage or imbibition (these modes 
20 respectively correspond to a saturation increase and decrease of the non-wetting phase). 
This phenomenon must therefore be taken into account to obtain representative relative 
permeability values. 
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The prior art concerning hysteresis effects in two-phase and three-phase media is 
described for example in the following publications : 

- Land C.S. : « Calculation of imbibition relative permeability for two and three- 
phase flow from rock properties », Trans AIME 1968, VoL243, 149, 

5 - Larsen J.A., Skauge A. : « Methodology for numerical simulation with cycle- 

dependent relative permeabiUty », SPEJ, June 1998, and 

- Carlson F.M. : « Simulation of relative permeability hysteresis to the non-wetting 
phase », SPE 10157, ATCE, San Antonio, Texas, 4-7 Oct. 1981. 

Figure 6 typically shows the course of the two-phase permeability curves 
10 resulting from drainage up to irreducible saturation in wetting fluid (M), then imbibition 
up to residual saturation in non-wetting fluid (NM). The hysteresis phenomenon occurs 
at two levels. At equal saturations Sg, different numerical values are obtained and the end 
point reached is an unknown that depends on the cusp point SgM from which the 
displacement mode is changed. This phenomenon is usually attributed to the trapped 
15 non-wetting fluid fraction. At equal saturations, we thus do not have the same quantity 
of mobile fluid, which distorts the flow characteristics. 

Practically aU the models taking account of hysteresis effects involve Land's semi- 
enq)irical relation : — = (1) 

where Cl represents Land's constant. This relation relates the initial saturation Sg, to the 
20 residual saturation Sgr in non- wetting fluid, in order to evaluate the saturations in trapped 
and free non-wetting fluid. Assuming that this relation is valid whatever the saturation, it 
is applied to determine the intermediate mobile fractions during displacement. In the case 
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of two-phase flows, associating this relation with permeability models provides 
satisfactory results. 

In the case of three-phase flows, the hysteresis of the relative permeabilities Krg takes 
on a particular form. A displacement hysteresis is e^erimentally observed, as it is the 
5 case with two-phase flows (extent of the direction of variation of the saturations), as well 
as a cycle hysteresis since the permeabilities depend on the saturation record. In Figure 7, 
the relative permeability curves Krg corresponding to a first drainage and imbibition cycle 
(Dl and II respectively) are distinct fi-om the corresponding curves (D2, 12) of a second 
cycle. 

10 The relative permeability model developped by Larsen et al takes these two forms of 

hysteresis into accovmt. Starting from an approach combining Stone's model in parallel 
with Land's formula and Carlson's interpolation method, an approach where only the 
displacement hysteresis is taken into account, Larsen et al have introduced an en^irical 
reduction factor that is a function of the water saturation, which allows to approximate 

15 to the permeability reduction of the gas associated with the cycle hysteresis. 

SUMMARY OF THE INVENTION 

The modelling method according to the invention allows to optimize faster and more 
realistically the displacement conditions, in a porous medivim wettable by a first fluid, of 
two-phase or three-phase mixtures including this first wetting fluid and at least a second, 
20 non- wetting fluid. It therefore provides operators with a more reliable tool for evaluating 
notably the best displacement modes of the fluids in the porous medium, in drainage and 
imbibition. It is based on a fractal representation of the porous medium with modelling of 
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the pores by a distribution of capillaries with a fractal section, considering a stratified 
distribution of the fluids in the pores, the wetting fluid being distributed in contact with 
the walls and around the second fluid (or around the other two in case of a three-phase 
mixture). 

5 The method comprises in combination : 

- experimental determination of the variation curve of the capillary pressure (Pc) in 
the pores of a sample of this porous medium in the presence of a wetting fluid and of at 
least one non- wetting fluid (by injection of mercury in a sample placed under vacuum for 
example); 

10 - determination, from this capillary pressure curve, of the fractal dimension values 

corresponding to a series of given values of the saturation in liquids; 

- modelling the hysteresis effects that modify the mobile saturations of the fluids 
effectively displaced as a fimction of the number of drainage and imbibition cycles 
undergone by the sanq)le, involving different non-wetting fluid trapping or imtrapping 

15 constants according to whether a drainage stage or an imbibition stage is carried out ; 

- modelling the relative permeabilities directly in the form of analytic expressions 
depending on the various fractal dimension values obtained, and 

- entering the relative permeabilities in a porous medium simulator and determining, 
by means of this simulator, the optimum displacement conditions of the fluids of the 

20 mixture in the porous medium. 

The method is applied for example for determining displacements of fluid nuxtures 
comprising a first wetting fluid, a second, non-wetting fluid and a gas, considering a 
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Stratified distribution of the fluids in the pores, the wetting fluid spreading out in contact 
with the walls, the gas occupying the centre of the pores and the second fluid being 
distributed in the form of an annvdar layer in contact with both the gas and the first fluid. 

The method can notably be applied for determining, by means of a reservoir 
5 simulator, the optimum characteristics of substances added to wetting fluid slugs injected 
in a formation alternately with gas slugs, in order to displace hydrocarbons in place, or 
those of a fluid injected into the soil in order to displace polluting substances. 

Modelling of phenomena by means of the present method has many advantages. It 
allows a better correspondence with the results obtained in the laboratory because the 
10 physical phenomena are better taken into account. The results of the model are therefore 
better in case of a scale change for example, for modelling an application in an operation 
field. 

The calculating time is reduced by comparison with the time required when tables 
are used as in the prior methods. Fractal type modelling can better deal with the 
15 hysteresis efiects encountered when using WAG type injection processes. 

The results of the method can fiarthermore be perfectly integrated into many 
reservoir simulators : 3D, heterogeneous, con^osition simulators, etc. 

E^loitation of the results by application softwares is fecilitated. It is no longer 
necessary to perform risky interpolations as it is the rule when working fi-om discrete 
20 values of the result tables in order to draw isoperms for exanq>le. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Other features and advantages of the method according to the invention will be clear 
from reading the description hereafter of a non limitative example, with reference to the 
accompanying drawings wherein : 

5 - Figure 1 illustrates, in the form of a table, the connections existing for a three- 

phase mixtxire between the experimentally obtained relative permeability values of a fluid 
and the saturations for two of the three fluids, 

- Figure 2 is a fractal representation of a pore, 

- Figure 3 diagrammatically shows the distribution of the phases of a three-phase 
10 fluid in a fractal pore with the wetting fluid W in contact with the wall, the gas phase G 

spread over the greatest part of the volume of the pore (radius Rq)* the oil O being a 
layer between the wetting fluid and the gas, 

- Figure 4 shows, as a fimction of the water saturation, the capillary pressure curve 
of a sanple of Vosges sandstone for example, whose local slope S is used to determine 

1 5 the fractal dimension of the pores, 

- Figure 5 shows, as a ftinction of time, the production curves obtained 
experimentally for gas (G), for water (W) and for oil (O), compared with the equivalent 
curves obtained by simulation by means of the method according to the invention, 

- Figures 6a, 6B show the effects of the displacement hysteresis affecting the 
20 relative permeabilities of the wetting flxiid and of the gas, and respectively, 

- Figure 7 shows the combined effects of the displacement and cycle hysteresis 
observed experimentally on the gas phase. 
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- Figure 8 shows the key cxirves used for modelling the untrapping constant, 

- Figure 9 shows the result obtained with the model within the scope of a WAG 
type injection, and 

- Figure 10 shows the validation of the model on a WAG type experiment for the 
5 recovery curves of the three phases. 

DETAILED DESCRIPTION OF THE METHOD 

The method according to the invention allows to determine the three-phase relative 
permeabilities of porous media by using a fractal type model of the porous mediiun, on 
the basis of an approach described by : 
10 - Kalaydjian, FJ-M et al. : "Three Phase Flow in Water-wet Porous Media : 

Determination of Gas-oil Relative Permeabilities imder Various Spreading conditions", 
68th Ann Tech. Conf. and Exh. of the SPE, Houston, Texas, 1993. 

The method according to the invention comprises, as mentioned above, modelling 
the flow of the phases with distribution of the phases in the fractal structure of each pore. 

15 In case of water and oil flowing in a water wet porous structure, the oil flows into the 
volume of the pore. In case of a three-phase flow, there is a stratified distribution, the 
water, which is the wetting fluid, flows along the walls of the pores, the gas circulates in 
the volume of the pore and the oil flows between the gas and the water. The saturations 
are calculated as the relative surface area in a cross section occupied by each of the 

20 fluids. 

At equilibrium, all the grooves having a radius greater than Rk, which is given by 
Laplace's law Pc = 2 y/R^, are occupied by the gas, and the smallest tubes by the two 
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other fluids (water and oil). The wetting fluid saturation is thus expressed as the fraction 
of the area of the occupied tubes. 

Calculation of the fraction of the area of the capillaries occupied by the water for all 
the radii between Rj^ and Roq leads to the following expression : 



and, as Pc = 2 y/R^, the correlation between the capillary pressure and the saturation of 
the wetting pheise is given by : 

1 

where Sw is the saturation of the wetting phase. 

10 The graphical representation of this correlation in a log-log diagram is a straight line 

starting from the point (S^,, Pc) corresponding to the largest capillary of the fractal 
structure with a radius Rq, One may suppose that : 

- radius Rq, which is first invaded when mercury is injected (Fig.4), corresponds to a 
saturation value of the order of l/r = 10^. Each segment of the capillary pressure curve is 
15 a part of a line starting from Rq (assumed to be the same for all the different segments), 
corresponding to the aforementioned correlation Pc, S^- Each line has a given slope, a 
fractal linear dimension can be associated therewith. The values of the slope range from 
-1.5 to -3.3 as shown in Fig.4, which leads to values of the fractal linear dimension Dl 
between 1.3 and 1.7; 

20 - each domain is reached by the mercury for saturations corresponding to the place 
where Rq is found on each line. 
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The saturation of the two liquids when the gas is present in the pore is calculated as 
explained above for a phase : 
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assuming that the two liquids occupy the tubes whose radii are less than or equal to Rj^ 
5 and the gas, the centre of each pore. The oil saturation is the relative area of the cross 
section of the grooves occupied by oil, whose radius is less than or equal to R^. 

Hysteresis modelling 

Calculation of the relative permeabilities requires determination, for each phase, of 
the fraction circulating therein and consequently systematic estimation of the saturations 
10 corresponding to the stagnant fractions. This must be done for the two cases studied, for 
example drainage of water and oil by gas, and water imbibition. 

The original feature of the modelling procedure proposed here is that it takes 
account of the hysteresis directly at its origin, i.e. at the level of the non-wetting phase 
trapping and untrapping phenomena. In Figure 8, if the curves could be deduced exactly 

15 by translation, this would mean that part of the gas has been trapped during the 
secondary drainage and does not take part in the flow. In fact, this is not the case but the 
fact that, at equal gas saturation, the permeability is lower during D2 than during II 
means that the trapped gas fraction is larger during drainage. In other words, the non- 
reversibility of the permeability curves can be explained by a non-reversibility (hysteresis) 

20 between the trapping phenomenon and the untrapping phenomenon. 
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Land's formula (Equation 1) is thus kept, but an untrapping constant, valid during 
the drainage stages, different from Land's constant valid during the imbibition stages, is 
introduced. It thus all comes down to modelling the evolution of the untrapping constant 
during the cycles. 

5 The characteristics to be taken into account in the formulation are as follows. We 

consider that Cp = Cl during all the trapping stages because we are directly in Land's 
conditions of application where a single constant is enough to describe the phenomenon. 
When Sgr is low at the end of the drainage process, Co is close to Cp. This simply means 
here that the gas is easier to imtrap when a small amoxmt thereof is present in the sample, 

10 because one can the consider that it is present in large pores that can be readily 
reconnected during re-injection. When Sg^ increases, Co becomes greater than Q, which 
shows a less eflScient untrapping process. 

For high gas saturations, a certain reversibility can be encountered, but a trapped gas 
fraction remains inacessible in small pores. A low single curve is thus reached, which 
15 corresponds to a pseudo two-phase case where the oil is no longer mobile. 

According to this representation, based on experimental observations and working 
hypotheses, the value of Cp passes through a maximum since the untrapping constant is 
equal to Cl when Sg^ is low and when it is maximum, when the low-mobility curve is 
described. The following expression allows to reconcile all the previous aspects ; 



Co- 



^ -K 

rg ^rg min 



Dl 



(2) 
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This formulation comprises several parameters : 

Cdm : it takes into account the difference between trapping and untrapping, 
Cl :Land's constant, 
K^mm " low-mobility curve, 
E : calibration parameter, 

K'^ : value of the relative permeability to gas at the beginning of the previous imbibition, 
K^^jg : value of the relative permeability to gas on the first drainage curve for the gas 
saturation corresponding to K^^, 

Whatever the order and the natvire of the cycle considered, relation (2) and Land's 
relation (1) allow to determine the trapped and mobile saturations by means of Land's 
formula while taking into account the three-phase character of the hysteresis 
(displacement and cycles) without using empirical reduction fectors. 



C is equal to Cl or Cd according to the displacement mode, 
Sgf : firee gas saturation, 
Sgj : trapped gas saturation. 

Calculation of the relative permeabilities 
A) Relative permeabilities of liquids 

Application of Poiseuille*s law to each capillary of the bundle for the phase that 
occupies it allows to calculate the water and oil relative permeabilities (Fig.5). 




andS^=S^+S, 



(4) 
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Other experimental studies (Larsen et al) have shown that there is a relation allowing 
to relate the residual oil saturation during the cycles to the trapped gas fraction. 

(Sor)sgt=o represents the saturation in residual oil left in place in the medium before gas is 
5 trapped. 

If we consider only the circulating fraction which contributes to the hydraulic 
conductivity, the relative permeabilities for water and oil are e3q)ressed as follows : 

R is the reduction factor linked with the trapping of the non- wetting phase. 

In these expressions, it is usefiil to mention that : 
10 - the irreducible water saturation Swi is assumed to be stable, 

- the size range of capillaries occupied by the mobile oil is calculated as the difference 
between the sizes of the capillaries occupied by the two liquids with the total liquid 
saturation Sl = Sq+S^ and those of the capillaries saturated in water and stagnant oil, 

- Kj.o(2ph.) is the value of the oil relative permeability determined by an imbibition test 
15 with water and oil. When only the water and oil phases are present and since the sample 

tested is water wet, the oil will flow through the section of the pore exactly as the gas in 
a three-phase flow. 
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B) Gas relative permeability 

Since gas is a non-wetting phase, it occupies the central space of the pore and it 
spreads towards the periphery thereof as the gas saturation increases, however without 
coming into contact with the solid wall (Fig.6). One considers that the gas circulates in a 
5 single pore whose radius Rg is given by the relation : 

Rg = R0+R1+R2+ +Rk. 
the gas permeability being then given by : 

K^-f^r,n^a-is,^s^ry (5) 



a = , Dt being the linear fractal dimension of the porous medium and Sl 

2-Dl 

10 the total liquid saturation equal to 1-Sg. 

The relative permeability model allowing to calculate Kj^, K^o and g is installed 
in a simulator such as ATHOS® or GENESYS©. This allows to caUbrate experiments 
carried out in the laboratory and to optimize the conditions to be satisfied in order to 
displace petroleum fluids in place in a reservoir, either by gas injection or by alternate 
15 injection of water and gas slugs (a method referred to as WAG), by taking account of the 
pressure and temperature conditions prevailing at the production depth. 

Validation 

The method according to the invention has been validated by means of various types 
of experiments : 
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gas was injected into porous media containing water and oil, under various 
conditions. It can be seen in Fig.5 for example that a very good accordance is obtained 
between the production curves of the three phases (water, oil, gas) obtained 
e^qjerimentally and those predicted by the reservoir simulator fed with the data obtained 
5 in accordance with the method, 

- gas and water were also alternately injected (WAG type injection). Fig. 10 shows 
that an excellent accordance is also obtained in this case for the three phases throughout 
the experiment. 



